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or Key^ 

Thompson and Stevens ji 3] proposed a model for 
emulsion polymerisation in a C.S.T.B., based on the popu- 
lation balance approach and instantaneous termination of 
radicals. Their model has been extended by integrating the 
■work of Sundberg and Min 1 22 I to calculate the mole- 

cular weight distribution, total number of particles, rate 
of polymerisation and the absolute particle size distribution. 
Chain transfer to monomer has also been included in the model. 
Effect of process variables viz; i) Initiator concentration, 
ii) Soap and iii) Residence time on these quantities has been 
studied. The theoretical predictions have been found to be 
similar to those of DeGraff 's |£oj and in some cases are 
better than the latter. Some experiments on the emulsion 
polymerisation of styrene in a C.S.T.P. have been done. These 
experimental results do not compare well with the theoretical 
predictions because of some inadequacies of the apparatus. 



CHAPTER 1 


INTRODUCTION 

Tracing the history of emulsion polymerisation 
process, one finds that it dates back to about fifty years. 
It has Reined great popularity in industry and has been 
to-date one of the most important processes for the manu- 
facture of rubbers, plastics, coatings and finishes.. 

The factors "which have contributed to emulsion 
polymerisation’s utility are i) easy operation of the 
reactor ii) good temperature control iii) low viscosity 
of polymer latex, and iv ) considerably higher molecular 
weights, at high reaction rates than that obtained by 
bulk techniques. 

The commercial techniques can be broadly classified 
into i) batch and ii) continuous systems. The batch system 
employs simple equipment and is extensively "used for low 
production rates and specialty applications. 

Batch processes have been exhaustively studied by 
Harkins [ 1 ] , Smith and Ewart [2] and Gardon [3] etc. The 
major theoretical and experimental aspects have been 
reviewed by Alexander and Naper [4] , Vander Hoff [ 5 ] and 
Min and Ray [ 6 ] . 
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The increasing demand of polymers has necessiated the 
use of continuous reactors and much attention has been given 
to their study. The scientific understanding of the conti- 
nuous process is still far from complete as it involves 
some conceptual problems that are absent in batch systems. 
Tor example, one such point of interest is the stability of 
continuous emulsion reactors. Gershberg and longfield [7] 
noted periodic behavior under steady inputs in their early 
experiments. Similarly Omi et al. [8] observed oscillatory 
behavior when polymerising styrene in emulsion. Satisfactory 
modelling of this behavior is yet to be done. Yet another 
problem is to find the equivalents of reaction regimes 1 , 

2 and 3 of batch systems, since in a well stirred system, 
micelles and monomer droplets can never disappear. 

In a continuous stirred system the various species 
spend unequal times in the reactor. The concept of resi- 
dence time distribution of the particles has been exploited 
by some workers. These investigators theorised that each 
particle is a small independent batch reaction unit and 
resides for different times in the reactor. A knowledge 
of batch behavior combined with the RTD is assumed to re- 
present the C.S.T.R. characteristics. The ma.jor difficul- 
ties with this approach are i) complex and long computations 
and ii) the difficulty of identifying the regime to be used 
for formulation. 



Hulburt and Katz's work [9] has generated interest 
in the application of the population balance technique to 
emulsion polymerisation, Funderburk [lO ] and Stevens and 
Funderbunk [11] illustrated the use of population balances 
in a G.S.T.R. These workers were followed by Thompson [13] , 
Stevens [12] ., Cauley and Giglio [13 ] . Their work will 
be reviewed in detail in the next chapter. 

Min and Ray proposed the most comprehensive popu- 
lation balance model taking into account a variety of 
physical effects but, did not generate numerical results [14] 
Ballal [15] used Min and Ray's model for his work. He 
could obtain number of particles, particle size distri- 
bution, rate of polymerisation and the molecular weight 
distribution but his model was too complicated and his 
calculations were time consuming. 

The present work has been done using Thompson and 
Stevens model [ 1 2 ] which is much simpler than Min and Ray's 
This model has been suitably modified to generate molecular 
weight distribution,, rate of polymerisation, the total 
number of particles and particle size distribution. A 
variety of physical effects have been accounted for. The 
theoretical predictions have been compared with the 
experimental data of BeGraff and Pohlein [16] . 
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Some experiments have been conducted on a C.S.T.R. 
using styrene as monomer, Potassium persulfate as initiator 
and Potassium oleate as the emulsifier. 
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CHAPTEP_2 

FHYSI_ CAI PROCESSES AMT MOD EL 

The constituents of the emulsion are i) water, 
ii) monomer, iii) soap and iv) initiator. Water forms the 
continuous phase. Monomer is generally insoluble in water. 

It is present mainly in the form of drops and to a minor 
extent in the dissolved state. Monomer also gets 'solubilised' 
in the soap micelles. The soap forms micelles and some goes 
to stabilise monomer drops and polymer particles. A neglegibly 
small amount goes into solution. The initiator chosen is 
water soluble and thermally dissociates into radicals. 

Some of the initiator radicals are captured by micelles 
and polymerisation is initiated. Thus polymer particles are 
created. Poe [18 ] proposed that the initiator may polymerise 
the dissolved monomer in the aqueous phase. Such radicals 
may interact to produce dead chains, which might precipitate. 
Such dead chains may coagulate and also form particles. The 
live radicals created in the polymer particles may desorb. 

Some workers proposed that live radicals desorbing from 
particles may react with radicals in the aqueous phase to 
precipitate dead chains. And as mentioned earlier these 
dead chains on coagulation may nucleate particles. 
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The initiator radicals can also enter polymer 
particles and generate growing polymer radicals-. They can 
also terminate existing live radicals in particles. Yet 
another phenomenon is the desorption of live radicals. Some 
workers have proposed that desorbed live radicals may re- 
enter the polymer particles. Due to the variety of physical 
processes there can be different number of radicals in 
different particles. Live radicals grow for different times. 
This generates a distribution in molecular weights. 

The monomer concentration in particles is maintained 
constant by diffusion from the monomer droplets while poly- 
merisation proceeds in the polymer particles. The live 
radicals can transfer their activity by reacting with a 
chain transfer agent or monomer to create a dead chain and 
a live radical. Due to the' growth of live radicals within 
the polymer particles, the latter grow in size. Since 
particles form at different times and grow intermittently 
there is a particle si.ze distribution. 

The monomer droplets and polymer particles are 
stabilised by adhering soap, (inspite of this there may be 
coagulation of the polymer particles). The micelles break up 
to provide the soap required for stabilising the growing 
particles . 
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A model then has to account for all these processes 
to predict the particle size distribution, number of par- 
ticles, rate of polymerisation and molecular weight 
distribution. All these processes have been incorporated 
in Min and Bay’s [14 ] population balance model. But they 
did not give numerical results. Most workers have generally 
neglected the following processes; i) particle coalescence, 
ii) chain transfer, iii) initiation in the aqueous phase, 
iv ) particle formation by precipitation of dead chains and 
coagulation in the aqueous phase and v) re-entry of live 
radicals into particles. 

Punderbunk [10] and Bunder bunk and Stevens [11] 
wrote population balances for a C.S.T.B. They calculated 
only the particle size distribution. They did not consider 
live radical desorption. They assumed that all particles 
contain the same number of radicals. They assumed the value 
to be that given by Stockmayer [t 0 )] which required the 
evaluation of an integral of a complicated ratio of Bessel 
functions . 

Thompson and Stevens [12] combined particle size 
variation with the Smith -Ewart [2'] recursion relation for 
radical balance in their population balance model. Gar-don’s [3 ] 
expression for particle growth rate was used and was assumed 
to be independant of particle size. A feed stream free of 
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polymer particles was assumed. Radical capture was assumed 
to "be proportional to the exposed surface area of the 
particles and micelles. They calculated the dimensionless 
number density of the particles as a function of their size. 
However the total number of particles, the rate of poly- 
merisation and the M.W.D. of the polymer product were not 
calculated . 

Thompson et al. [ 1 3 ] used the model of Thompson and 
Stevens [ 12 ] (TS Model). They assumed instantaneous termi- 
nation of live radicals in particles. This implies the 
presence of either zero or one radical in the polymer 
particle. They could obtain analytical solution of the two 
population balance equations framed for the two types of 
particles. The}? - calculated the non-dimensional particle size 
distribution and compared these with Funderburk 's experi- 
mental data [W>M] . However the criticism levelled against 
Thompson and Stevens work remains valid here too. The 
absolute value of the total number of particles could not 
be calculated since the variables in the balance equations 
were made dimensionless using -NT the number of zero-size 
micelles, which could not be estimated. For the same 
reason rate of polymerisation could not be calculated. 

It is to be noted that T.S. model is incapable 
of computing M.W.D. 



Recently Sundberg [17 ] in his hatch study of emul- 
sion polymerisation used population balances. He made var 
ables dimensionless by choosing measurable quantities. 
Along'witb the population balances he wrote a material 
balance on soap and solved a set of simultaneous partial 
integro-diff erential equations numerically to obtain the 
absolute particle size distribution. He also obtained the 
polymerisation rate and conversion as a function of time. 
Radical desorption was not considered in his model. 

The present model combines the T.S. model with 
Sundberg's model and Min’s [ 22 ]work for the study 
of emulsion polymerisation in a C.S.T.E. Particle size 
distribution, number of particles, rate of polymerisation 
and M.vV.D. have been calculated. 
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CHAPTER 5 
MODELLING EQUATIONS 

The following assumptions have been made in this 

chapter 

(i) It has been assumed that as soon as a primary radical 
enters the polymer particle, instantaneous termination of 
growing radicals occurs. This implies that only two types 
of particles are present in the emulsion; a) particles 
which contain one live radical and b) particles which do 
not contain any radical. Radicals and dead chains can be 
of various length and are housed in the particles. 

(ii) The micelles and polymer particles are presumed to 
be the sole sites for all reactions except for the disso- 
ciation of initiator molecules. Capture of initiator 
radicals is assumed to be proportional to the surface 
area of micelles and polymer particles. 

(iii) The reactor feed is assumed to be free of polymer 
particles . 

(iv) Particle growth rate is assumed to be independent of 
particle size. Gardon's expression; = Ki 

is used where K is the growth rate of a particle con- 
taining a single radical, i is the number of radicals per 
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particle and v is the volume of the particle. 

(v ) The initiator radical has been treated as a polymer 
radical of unit length. 

(vi) All radicals in the particle are assumed to react 
with equal reactivity independent of their length. 

The following kinetic scheme has been adopted 
K d o 

^2 y '21* (Dissociation) 

K T 

I* + M — £.j (Initiation) 

K 

? n + M P^ (Prorogation) 

K 

+ M _ T r . M + P^ (Chain transfer) 

K, 

P +1* , ^ M (Termination) 

n n 

Population Balance Equations Po r Determining Size Distri- 

bution, N m and R 
T _g 

let f(i,v)dv be the number of particles containing 
i radicals lying in the volume range of v to v+dv per cc. 
of emulsion. Note that i = 0 or 1 only, since instan- 
taneous termination has been assumed. Note that f(0,v m ) 
does not denote micelles. 
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Following Thompson and Stevens [12] the population 
balance equations can be written as:- 


2 


t ^ X f(0,v)] = 0 


= R ia a p N a [ f(l,v)-f(0,v) ] 




f (i ,v ) - -- ^ 27 ) 


L.H.S. is zero because ^ is zero for particles which do 
not contain any radical. The first term on F'.H.S. represents 
a gain and a loss term through radical capture. The second 
term represents a gain through desorption of live radical 
and the last term shows loss through outflow. Using 
dimensionless variables and parameters defined at the end 
of the chapter, we have:- 


X 

o 


(A* y- 1/3 ♦ A, y^ 3 ! 


fl 

A y 2/3 7 - 


( 1 ) 


Similarly 

[Kf(l,v)] = R ia a p \ [f (0,v ) - f(l,v)] 
K a 

- — ~~ E f(l,v) - f (i .v ) 
v x 

Non-dimens ionalising 
dX 

-3yL = a 4 y 2/3 (X 0 - x i> -^i x l - A 5 X 1 

173 

•tr • 


( 2 ) 
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Adding (l) and (2):- 


dX 1 

w 


- VW 


(3) 


Solution of this equation yields; 


X 1 (y) . 

ln ZjTTT = 2 A 5 (l -y ) ~ a 5 v > 


3xAgA 


2 A, 


5. 


In 


1+A 3 y 

1+A, 


2/3 


' ^'[tan^UjI/ 2 y 1/3 ) - tan' 1 (a, 1/2 )] (4) 

A 3 

where A^ = A^/A^ 

A b.c. has to be written to find X^(l). for th.is, the 
procedure given by Sundberg [17 ] is followed. Writing a 
balance in a volume of width ev m (where e is a small 
number), we have : - 


x f(l,v m +ev m ) = a p U A t f (0,v m + dls)-f(i , Tm+ -J!)] e v 


-pr- 

m 2 


m 


K a £ v 

~ 2 f(l ,v + ) e V 

v m 2 m 


f (1 ,v + fV ) ev 


+ R . m a N . ( 5' 
la -urn A v 


The growth rate terms have been written for particles coming 
in and leaving a volume element of 6 v ffi at v m and v m + e v m ’ 
The first, second and third terms have to be written for 
all particles lying in the volume range of v m to v ffi + e v m * 
Hence their volume has been taken as the mean valpe in this 
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range, i.e. v ffi ev ffi . Here only one growth term (output) 

+ 

is present as a particle of volume v„, can not exist. The 

m 

last term appears because of creation of particles from 

micelles at volume v . 

m 

Taking limit as evr^O 

m 

Kf ( W = ^ia m a m 

or X 1 (1) = h -f- (6) 

0 

Material Balance For Soap: In a well stirred C.S.T.E., the 
concentration of species in the reactor is the same as the 
concentration of the species in the effluent stream. The 
polymerisation rate is proportional to the concentration 
of X 1 - type of particles (see eq. 13). The concentration of 
X-j-tjrpe of particles is proportional to m (from 6). Hence 
m, the number of micelles in the out flowing stream can 
never be zero, for the model to predict non-zero reaction 
rates. As long as the micelles are present, the soap 
concentration in the aqueous phase will be equal to the 
saturation concentration in water. Hence a material balance 
on soap would mean equating the soap supplied through- 
micelles to the reactor, to the soap leaving the reactor in 
the fnrm-of 'soap' on particle surface and monomer droplets, 
and, the free micelles in the output stream. 
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2/3 2/3 1/3 


m o a m = m a m + A p + v d 3 } 


(7) 


Monomer droplets though larger in size than micelles, are 
far fewer in number. Thus soap required to stabilize them 
may be neglected. 


or m a = m a 
o m m 


f * r 

+ Jf(o,v) a p d^+ /f(l,v)a dy 

J r 


on simplifying; 


1 = E + 

in 


cc 

r 2/3 
x 0 y dy + 




2/3 


x -,y dy 


Prom (3) we haver 

(X 0+ I )y 2/ by = - l 


w 

' dX 


ay y 2/3 fl y 


( 8 ) 


Simplifying after using X^ ( ) = 0, 

m 1 


1+ ii + ^TW3^ 

A c 3 A 


L 5 J 




X 1 

Substituting for “(-]) (4-), 


m 

m. 


(9) 


can be calculated 


oo 


Now N^ = 


[f(o,v)+f(l,v)J dv 


(10) 


'Vi 05 


= m. 


(X Q +X 1 ) dy 


Prom (3); N n 


m 


po 


dX, 


dy 


dy 


A 4 m 

Using (6), we have; N = — x -- x m 

5 o 


(in 
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Thus N t can he calculated. 


Rate of polymerisation, R is given by;- 

P 

K (M ) x radical cone. 

r = — E. — E __ 

P N a 


But radical concentration is given by the total number of 


particles containing one live radical/cc. of emulsion. 
K„(MJ /°° 

f ( 1 ,v)dv (12) 


R_ = x 


P 


N 


A 


W A-m 

X “m” X m 0 


A 


o 


o 

X 


1 


xgry ay 


(13) 


Therefore R^ can also be computed. 


Particle size distribution can be obtained as follows; 

The cumulative size distribution is given by 

[ jf(0,v)+f(l ,v )] dv 
P(y) = ^ 


|(0,v)+f(l ,v)j dv 


(14) 


Using (3) and simplifying: 

, , x i^y) 

F(y) = 1- y 


(15) 


Thus particle size distribution can be obtained. 


Population Balance For MWL Determination 

Let L(n,v )dv be the number of particles lying in 
the volume range of v to v+dv , having a radical of length n. 
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>T [KX.(o.t)] = - E ia a p Nj I(n,v) - L (n,v) 

+ E ia a p % f(o ^> & n1 +K p(« p ) [Kn-1,v)-L(n, T )] 

~ My) _ K tj ,(M p )i(n,v) + K tr (K p ) f(l,v) 6 p1 (16) 

The fourth term on R.H.S. represents a gain and a loss term 
through propogation. The last two terms have been written 
to account for chain transfer to monomer. 


N on-d imens i ona lis ing (16), we have:- 


f y *'Cn,y> = - l’ - 1 % A 4 y^ X Q 6 


-1/3 


2/3 


nl 


+ A 7 [ l*.(n-1 ,y)-L* ] A q L*+ A g X 1 6. 


nl 


(17) 


Following a similar approach as in (5) and writing of 
balance on a volume range of v m to v m + &v for l(n,v); 

♦ 

0 = Q1 (n,v_ + e? m ) £v m . + KL (n,v + ev ) . 

m m m m m R 

2 

- V B f(0,T m + ) ev m E. a a p H A » n1 - V R B ia H A a m m 6 p1 

2 

- K tr (M p ) 1 ( n ,v m + ) ev m . V R + K t r (M p )f(l ’ v m + ) x 


ev m'Vn1 + V R L (n ’V fjm 5 ‘ + r^ 3p ^ ) ( 18 ) 

l.H.S. is zero because there are no live radicals in the feed 
and ^( n > v m ) is zero for the same reason given for f(l,v m ). 
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Talcing limit as ev m --f 0 


K L(n,v ) = F m a N. 6 

m ia m A nl 


or 


and 


K Un,v m ) = Kf (1 ,v m ) 6 n1 
L(1, v )= f(l,v„) 


m 


l(n,v m ) = 0 for n> 1 


' . L (1,1) = X,(1 ) 


(19) 


Let G^(n,v)dv be the number of dead polymer chains of 
length Inin particles containing i radicals and lying in 
the volume range of v to v+ dv . Minimum value of n for G^(n,v) 
is 1 (obtained by the chain transfer to monomer by an ini- 
tiator radical) and for G Q (n,v) is 2 (obtained by termination 

of an initiator radical). 

G (n,v) . 

-S-j G 0 (n,v) R ia a p U A + [ «n-1,y) 


0 = 


+ G^(n,v) ] R ia N^ + G^(n,v) K Q 


( 20 ) 


Non-dimens iona Using (20); 


0 = - A 5 G* - A a y 2/3 (G* - G*) + A^G ^ + A 4 y 2/3 L (n-1,y) (21 ) 


173 


and , 
... 2 __ 


G t (n,v) . 

[KS^n.v)] = - T + G o (n ’ T) E ia a p N A+ K tr (^)X. 


*“** Gj R . -fet Na 
1 ia p A 


G 1 K o a p 


v 


( 22 ) 
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Non-dimens ionalising (22):- 


dG 


_1 

cy 


2/3 


* * 


- A 5 ^ + -G 1 )- A 6 G 1 + A 8 I (n,y) 


(23) 


y 


1/3 


Writing a balance en G (n,v) in the volume range of v ffl to 

v m + ®V 

0 = G o (n ’ v m + !lin } " F: ia a p N ;J L(n_1 ’ v m + )+G 1 (n}V m + flir ; 

2 T 2 2 

+ R,.„ a„ U, a v„ G„(n,v„+ ev m ) - ev, Gl (n,v m + ev m ) 


m la p ±i m o ’ m m 

2 


(24) 


OP i* i. 0 

v m ^ 

The second term on R.H.S. represents a gain, first 
through radical capture by l(n-1 , v ffi + £v m ), which creates 

a dead chain of length n and then by radical capture by 
G 1^ n ’ v m + ^ creates a G Q type of chain of length n. 


Talcing limit of £v 0 


m 


K o a m ® 1 (n ’ 7 m ) . E ia a m % [I(n '1 > V + G t s o ( n > T in ) 

“ + ^ 


V 


+ 3R. a N. Q (n , y ) 
la m A o m 


(25) 


On 


taking a similar balance on G-j(n,v) we have:- 


0 = G 1 (n,v m +ev m ) a v m + K G, (n,v m + a v m *t G, (n,v m + ^ ) a v 


+ E ia a p 


H a G^ (n,v m + % ) % - E la a p N a G 0 (n,v m + j-Vj,) e v m 


+ K tr (M >l(i,y « m ) ev m 

2 


« O C «l* ' .'.(.26 ) 
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Taking limit, ev^ =, o, we have:- 

S ( n,v m = ) " = ^ ^ or an y n * 

P* S N f* ( 1 T7 1 

Now from (25) :- G (2 ,v ) = — A§ — S — 4_ LiL — 

4 + E. a N. 

* ia m ii 

and G 0 (n,v m ) = 0 for n>2. 


Solution of (2?) and (23) would yield the values of 
* * 

G q and G.| for various values of y. Hence a plot of 
[*° 150 

I n(G Q + G^ ) dv / pn(G 0 + G^ )dv vs . n could be made, which 

V rrJ VfC 

would constitute the M.N.D. of the polymer product. 



m 


Jr-. 


'P 


m 


(M ) 
P 


“d 


b 

b 

m 

n 

R 


la 


v 

v 

i 

v 

c 

V 

9, 

x 

6 . 


m 


E 


m 


ig 


N 0MENCIATT3HE 

surface area of each micelle, cm 2 

surface area of particle of size v, cm 2 

total surface area of particles, cm 2 /cm^ emulsion 

density of polymer, gms/cm 5 

density of monomer, gms/cm^ 

number of radicals in a particle, dimensionless 
concentration of monomer in the particles, moles/cm/ 
rate constant for dissociation of initiator in 
aqueous phase, 1/cms 

rz 

rat<= constant for propagation, cm /mole .sec . 

rate constant for radical desorption, cms/sec 

number of micelles, ■' /cm'', emulsion, i^ub- 

chain length of polymer, dimensionless 

rate of absorption of radicals in particles, and 

2 

micelles, moles/sec. cm 
time, sec. 

volume of particle, em^ 
volume of micelle, cnr 5 

3 

volume of monomer droplets, cm 

3 

volume of reactor, cm 

monomer volume fraction in the particles, dimensionless 
mean residence time, sec. 

Kronecker delta function, 1 for i=g, 0 for i/g 

1.1 i *’■ -f-PUH 

mm 





It**. 



f(i,v)dv number of particles containing i radicals and 

lying in the volume range of v to v+dv/cnr emulsion 
L ( n , v , ) dv number of particles containing a live radical of 

length ,'n and lying in the volume range of v to v+dv 
per cc. emulsion 

G i (n,v)dv number of dead chains of length lying in particles 
of volume ranging from v to v+dv and containing 
i radicals/ cc. emulsion 

rate constant for chain transferwithmonomer 
moles of soap/cc. solution 

2 

surface area of soap per mole, cm /mole 
number of micellse/cc. of inflowing solution 


K 


tr 


tti. 


y 


u 


A 


Y , dimensionless 

.VagadroS 
•p (i,v) , Vm/^o 
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C HAPTER 4 
METHOD OF SOLUTION 
determination of and 

From ( 9 ), m/m can be calculated once the integral 

'’■'r-j'*, 

in the denominator I /X . ( 1 ) x y^^dy has been evaluated. 

'V '( 

For E , \X./X 1 (1 ) dy has to be evaluated. Both those 

ir ’ * 

integrals were evaluated using a 5 point Legendre integration 
scheme, after converting the limits to the range of -1 to +1 , 
using the following conversion: - 

fc r i 

if(x)dx = s jf(n)dn where n = — 

> : v 

Here a = 1 and b was an arbitrarily chosen large number. 
Convergence was tested by using different values of b till 
the value of the integral remained unchanged. See appendix II 
for the program used to integrate the expression. 

Mol e cular Weight determination 

As it would be tedious to calculate all individual 
values of L(n,v) and G i (n,v), Minsand Hay's proceedure 
of using moments for generating M.V.D. has been followed 
here. The moments of the L distribution along with the 
moments of the G^(n,v) distribution has been used to generate 
the moments of M.W.L. These moments can then be used to 
generate the M.W.L. 
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Let L be the z-transform of L (n,y). 


. ' . L = j> Jj (n,y ) z 


.n 


Multiplying ( 1 7 ) by z n and summing from n= 1 tooO we have. 


^ = -1 ( A^+Aq+A 4 y 2//3 + A^y "^ // ^) + z(AgX^+ A^X 0 y A^Lj A^L 

1 ,the n^ 'moment of L*(n,y), can be obtained from the 
z-transform, as follows:- 


u = <« fp n i 


n 


z=1 


The differential eauations for the moments of L(n,v) 
distribution can then be written as:- 


dl. 


dv “ fl n + liC o 1 o A 7 + ^l 1 ! A 7 + C 2 A 7 1 2 + 


+ nC n-1 A 7 1 n-1 +e 


where f = ~(A 8 +A 8 +A 4 y 2//3 + A g y 1//3 ) and g = AqX 1 +A 4 X Q y 2//3 
The b.c. for these equations are as follows :- 


1 . 


1 1= 1 2 = 1 3 = .... = l n = X^l), at y 


( 30 ) 


It is possible to obtain the moments of M.W.D. directly 
without computing the moments of G^(n,v). 

Adding (21 ) and (23) we have:- 


dG 

¥ 


-rr 

1 = -A 5 (G* + G* ) + A 4 y 2/3 L*(n-1 , y)+A Q L*(n,y) 


( 31 ) 
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«C 


Let 


A 


.n r (G + G 1 )dy 


th 

Then I, is proportional to ti e k moment of 


l*i © - • 


Then from (31 ) * — L^. 


i k j 
; sr n j 

Ar 1 


dG* f 

dy~ + 


1 , L 2/3-* 


* *, j 

+ AgL (n,y); f dy 


A 4 y (n-1 ,y ) 

_.(32) 


¥e have seen that G^(l) = G-. ( ^ )= 0 


or 


or 


/ iL | 

2 \ o /■* * x * i 

I k = J^-A 5 j A 4 y^° L (n-1 ,y ) + AgL I dy 


Sc 


1 i 

A 5 ! 


A.y 


(m+1 ) k L*dy+ ^8 j //n k L dy 




A 


c; -i 


” A = V A s j> 3 *^ k! 


p- Cic-p); 


m 


p L*(m,y)dy + ^8 j l.dy 


An 


Defining K v as=- 

& , 


K 


'k 


Ar 


j'~2/3'> < ^ k ; m p L (m,y )dy + Ag/Ar jl k dy 

> y k_ p F (k-pT • ' 


"b.c. K k = 0 , at X = 1 and K k = I k 

lim x vJC> 

K, on differentiating twice w.r.t. x yields a set of second 
K 

order differential equations which can "be solved using the 
eauations (29) and (30). This would "be a simpler method of 
obtaining I k than "by integrating (32). 
dK, A , 2/3 £ k! • A 


dx "A* fS pf Tk-p7j 1 P 

J f ' O 


l-_ (x ) + T“ li,(x) 


a 5 x k 


(33) 
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dx 


2 

„ d jK 

& • 

- 2 
dx 


At x 2/3 1(l) /ii ( x ) 

A x 1 o ( ' x; + A 5 ± o^ x; 
dK 


(34) 


# 

dx 


_ f+ Vt^rrti~J 


x 


2 _i/3 ! , il 4 _s ’ ^*8 \ 

f X 1/3 j + g(— - A - — ; 


,2/3. 

X 


.. .(35) 


dK 1 _ H 2/3 
dx i'- 5 


( l o + V + I ^ 1 ! 


& d 2 K. 


dK 


dx 


GX 


1 / 2-1 2-1 
( ^ x ~ 5 x x 


A x 2/3 A 

2 4 X _3 . „1 

,2 + C o 


3x 


(a 4 x 2/3 + a 8 y 


+ g 


V ' dK o 

T 2jT’ +f ^ + 

V +A 8 

2A 4 x 2/3 + Ag 


dx 


A, 


(36 ) 


(37) 


1 

■where C n = Ay 


dK 


1 = li x 2/3( 1 0+ 21 1+ 1 2 )+ A s 1 2 /A 5 


dx A 


(38) 


, d 2 K dK 9 2t dK. a 1/3 + 2 

& --4 = — ~{f+ 3i ) + aX (- 4y " ~ ^ 


dx 

5 

dx 


dx 


2 A1;y - ^/3 

+ c 2 + ) 


+ ffa ( 4|£ _ 2t 2 x 1/3 + 2tx -1/3 /3 + 2tc 1 x 2/3 . c 2 x 2/3 t+C 2 ) 


+ if ( ’a 8+ «4 x2/3) 

9 


* 0*0 


(39) 


dK<? A a o /•? x ^3 

ta 5 - 4 - 2/3 (l. + 31 1+ 31 2+ l 3 ) + jf 1, 


( 40 ) 
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& aS 

ai 


dK 


dx 


( 3 C 4 x 2 / 3 t + 3C 2 x 2 / 3 t + C 3 - C 3 t x 2 ^ 3 


+ 2t 2 C 3 x^- 3C 2 t 2 x 4 ' /3 - 0 3 t 2 x 2 / 3 + 4xt 3 + 


3. 2 x 


•1/3 


± - 4t 2 x 1/3 


- 4x 5 /V + 4xt 3 - ) + |/C C 3 + 4xt 3 - 6t 2 x 1 / 3 + 2x - 1/3 


3 tx 2 / 3 + 3G?x 2 / 3 t) + ^2( c 3 + 2 tx- l/3 - 2 t 2 x 1 / 3 ) + ^ (f+ 


20 ; 


tx’ 


'1 


+ ( At S + 7g 

V *5 


d^ 2 

A>2/3 

) 


cx 
(41 ) 


2/3 


AA = 1 ( V Jllfl) + 41, A x 2 / 3 + 61, x 2 / 3 + 41 £ x 2 / 3 

dx 4 A,- 3 Ac 2 iir ^5 


+ £ * 2/3 4 

5 


(42) 


d^K 

& 4 


( f +2 V -1/3) + -L (4A V + 4 H4 fx 2/3 
dx 2~- "4 f -T^t + 3 A 5 X j + 1 3 UA 7 A 5 t + 4 A 5 1X 


8 -1/3 


6 A rv Ai , 


A 


) + 1,( -xp + 12 A X 2 / 3 4 7 + 4x^ 3 ^ 


+ 3 x t: j " X 2 V _ Act ‘ A 


+ 6f ti x 2/3) + ( 1VA + x 2/3 + x 2/3 

Ac 1 Act A 5 5 


+ 3 x 


8 ,-1/3 _4 + x 

A 5 A 5 


4 A 4 2/3 „n , , / Vi , ViV „2/3 


+ 


Ar 


X 


+ 1 o ( Act 

+ JAA x 2/3 + X 2/3 A. + ^ fx 2/3 + f x" 1/3 

Ac Ac I ^ J "'5 


6h,A 


gA 4 15A 4 2/3 
A 5 t 4 A 5 X 


( 43 ) 
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wh ere 


t 


2A 

" + Ag 



C 






The h.c. are as follows:- 


fa 

dx 


V A 8 x (1) 

A 5 x i u; 


■when x = 1 


dK 1 

dx~ 

dK 2 

dx 

dK 3 

dx 


2 V Ag ^ 

Ar f 

5 


ihlis (1) 
t *,111 


8A a + An 

1,(1) 


when x = 1 


when x = 1 


when x = 1 


dK 4 

dx 


1 6A , + An 

1 1 X. (l) when x = 1 


Equations (34) to (43) have been solved by using a fourth 
ord ero Gill’s modified technique of Runge-Kutta method of 
numerical integration. 
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CHAPTER s 
E XPERIMENTAL ¥ QRK 

Ajg£ara_tus . - A cylindrival glass vessel 'with a hemispherical 
bottom was used as the reactor. Three stainless steel baffles 
were used. The reactor had a side outlet which allowed a 
volume of 500 cc.s. of the emulsion to be present in the 
reactor at any time. A stainless steel lid rested closely on 
the flange of the vessel as shown in the figure 8. An adopter 
was fitted centrally on it through which the stirrer rod was 
introduced. This was done to present leakage of the Nitrogen 
blanket and other vapors. A Remi stirrer (l/S E.P.) was used 
to stir the reactor contents. Two stainless steel tubes 
passed vertically down through the lid and were used to feed 
the reactor with the monomer and soap solution respectively. 
The feed rates were adjusted using needle valves. The tubes 
were connected to two reservoirs filled with soap and monomer 
respectively. The levels of both in the reservoir were main- 
tained to be constant manually. The reactor contents were 
heated by heating coil wrapped around the base of the reactor. 
The reactor rested on a screw jack so that it could be raised 
or lowered as desired to maintain proper head required for 
flow. The reactor contents were kept in an inert atmosphere 
by using nitrogen gas at one atmosphere pressure. 
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Experimental Procedure : 

The styrene monomer was washed repeatedly using a 
15?'° solution of ft. OH., to remove the inhibitor present in it. 
This was followed by washings with water to remove the last 
traces of ft OK from the monomer. Next the required weight of 
the initiator (potassium persulfate) was dissolved in the 
soap (potassium oleate) solution. 

The soap solution and the monomer were then allowed 
to flow into the reactor at the required flow rates. After 
some time, stirring was started. When the level of the 
emulsion in the reactor stabilized, heating was started. 

Once the temperature attained the steady value of 70°C, 
samples of the exit stream were collected. This was done 
for three different flow rates of the reactor feed such that 
the residence times were i) 25 mins., ii) 40 mins, and 
iii ) 50 mins . 

Polystyrene was precipitated from the emulsion using 
methanol solution. It was then drl®d in a dessicator and 
then was analysed in a gel permeation chromatograph using 
tetrahydrofuran as the solvent. 

Limitations of the Apuara_tus ? - 

The major limitation of the set-up was the stability of 
flow through tbp reactor. The outflow courd not be kept 



absolutely constant for a period of time . Hence the outflow 
from the reactor could not he maintained equal to the 
inflow to the reactor. 
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Oxidation of the monomer could not he avoided totally, 
hence there were some oxidation products like aldehydes and 
ketones in the reactor feed. 
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CHAPTER 6 

RESULTS A nd DISCUSSION 

Nu mer ica l Derails • The numerical values of parameters 
reauired for the evaluation of and A Q were taken 

from De Graff and Pohlein £ 1 sj , Min and Ray [£J and 
Brandrup and Immergutj^ 2.1 J . A list of all such values is given 
in Table 1 . Values of Ag were assigned arbitrarily as 0.05 and 
0.5. The cumulative particle size distribution, the total 
number of particles, the rate of polymerisation and the molecular- 
weight distribution were calculated as a function of the following 
process variables : i) Initiator concentration, Cl), ii) Sur- 
factant concentration (S) and iii) residence time (-). 
Calculations were made only for those condition as were reported 
by DeGraff and Pohlein 16 in their experimental study to 
facilitate comparison. 

Numerical computations were performed on Dec-10 computer. 
Determination of N the total number of particles, and Rp, 
the polymerisation rate required a minimum of about lOsec.s 
to a maximum of 1.2 mins. The evaluation of the moments of the 
molecular weight distribution, MUD, required a maximum of about 
2*5 mins, of C.P.U. time. 






31 


P arti cl e Size Bistributi nn 

No desorption* I 1 , the cumulative size distribution 
function has been plotted against particle diameter, D for 
different residence times in fig.1. The plot shows that as 
residence time is decreased, the distribution shifts towards 
smaller sized particles. This is to be expected since as 
residence time decreases, the particles grow for shorter 
periods . 

The effect of desorption of live radicals is to dec- 
rease the particle size. When desorption occurs the effect 
is to lower the time of growth of a particle. Hence the 
distribution is shifted towards smaller siz'ed particles. 

When desorption is operative. 

A good fit of the theoretical curves with the experi- 
mental data is observed at larger values of D, whereas some 
deviation is observed in the smaller size region. The theory 
predicts a smaller number of smaller particles than observed 
in the experiments. The present theoretical prediction of 
the size distribution are very close to that of DeGraff’s 
for large residence times. But the present theory shows a 
better match with experimental values at low residence times 
The effect of desorption is not appreciable for Ag= .05. 
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Since the effect of desorption is to decrease the particle 
si7.e, a slightly larger value of Ag may fit the experimental 
data better. 

and R p 

Effect of Initiator Concentration: It is observed from fig. 2, 
that as (I) goes up from about 0.3 (gm/100 gms.HgQ), both 
h T rp and Rp remain unchanged. The constant value reached for 

1 R 

N^, is about 2.14 x 10 " = /cc., and that for is about 
8.8 moles/lit , hr . It can be inferred that (I) does not 
affect Ny because after a certain threshold concentration of 
(I), the radical concentration in the reactor is fairly high 
in proportion to the micelles present and therefore no new 
particles are formed. Only the frequency of radical capture 
is increased. too is proportional to the concentration of 
growing particles. Thus R^ also reaches a constant value 
beyond a certain value of (I). However very low (I) would 
decrease both and B^ since then radical concentration 
would be very low. The theory also predicts that at zero 
concentration of (I) both and R^ become zero, as expected. 

When desorption is operative the theoretical curve for 
shows a gradual decline from a high value, till it merges 
with the line for no desorption. R shows an increase from a 

if 
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low value to the constant value attained by R fox no 
desorption, 'hen (I) is large, the radical concentration 
is high and hence the effects of desorption become unimportant. 
Thus at large values of (i) the curves for no desorption and 
with desorption should merge. The behaviour at lower values 
of (i) can be explained by noting that desorption impedes 
the growth of polymer particles. The smaller particles require 
less soap for stabilisation, hence more soap is available as 
micelles. This causes an increase in the number of particles. 
But as (I) increases and more radicals are available, they 
begin to counter the effect of desorption by promoting growth 
rate. Therefore Nrp starts decreasing with (I) and finally 
merges with the curve for no desorption. 


R is proportional to the number of growing particles 

i e X 1 H T , x i 

* -V-— t" v— .Now from 1 , can be written as equal to 

x o + x 0 + x 1 

A 5 +A 4 y 2/3 

— — —— - y . This function has been sketched -with y 

A 5 +2A 4 y Z/i +A 6 y 1/5 



y 
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It is seen that increasing desorption causes a corresponding 
X 1 


decrease in 


V x i 


for all sizes of particles. Therefore 


F’p decreases when desorption is present inspite of the 
enhanced values of N^, . It increases continuously with (I) till 
it begins to merge with the results for no desorption. 


Good fit of theoretical predictions for no desorption 
with experimental results is seen. It can he concluded that 
desorption is either inoperative or present only to a minor 
extent. 


Effect of Pesidence timet It is seen in fig. 3 that 
both and fall continuously with C . As Z is increased 
the particle size increases and more soap is required to 
stabilise them. Consequently less number of micelles are 
available for nucleation of particles thereby decreasing 
N^. Thus and R , which is proportional to N^, both 
decrease with ^ . However the theory predicts that when 

is equal to zero, both and would be zero. Thus both 
and R must first rise to a neximum value and then begin to 

Jr 

fall with "~C • This is to be expected because at low residence 

times, as the residence time increases more particles would 
be created and hence both and Rp would increase. 
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As discussed earlier desorption increases the number 
of particles and essentially leaves R^ unaffected* Thus the 
nature of variation of N^ and R^ versus remains -un- 
changed due to desorption. The R^ curve for finite desorption 
lies a little below the one for no desorption. This also is 

Xi 

to be expected as P is influenced by w- ,- v which as shown 

P 

earlier is responsible for lowering the value of Ft * 

The theoretical predictions without desorption show 
some deviation in the high T region, for both N n and R . 

JL .P 

The fit with experimental data is reasonably good for low- 
residence times. Some desorption effects can be detected in 
the plot. 

Effect of surfactant concentration: R is seen to 

Jr 

rise sbeeply with S in fig. 4- This is understandable because 
as more soap is fed more particles are created through 
micelles. Therefore both N^ and R^ increase with (Si. 

A fair fit with experimental data is found and the 
match is no better than that of DeGraff ' s theory with the 
experimental points. 

Number average molecular weight 

Effect of (I):- is observed to fall steeply with (I) in 
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fig. 5 after a certain value of (i) (around /gm /100 gms 
The rate of decrease is a little slower below this concentration. 
An explanation for this phenomenon is that, is zero when 
(i) is zero hut would increase with (I) for low values of 
(i). But as (i) increases beyond a certain value, rate of 
radical supoly to the particles becomes high and hence live 
radicals created within the particles are terminated faster, 
with the result that time available for their growth is less. 

Then small chains are produced and a low M n results. 

When desorption is present the curve for Ag equal to 
.05 follow the line for no desorption closely. However, the 
curve for Ag = .5 clearly show a maximum but follows the 
curve for no desorption at higher values of (I). Both the 
lives (ie for Ag=0.05 and 0.5) are below in for no desorp- 
tion. This can be explained by the fact that at high (i) the 
desorption effects are not iimportant since radical supply 
is large and both curves with and without desorption should 
follow each other closely. Again at very low values of (I), 
the radicals enter the particles less frequently and desorp- 
tion drastically decreases the growth time of the radicals 
leading to a low value of M . 

The fit of the theoretical prediction with experi- 
mental data is reasonably good for Ag= 0.0 and Ag=0.05. 
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It can be inferred that desorption effects are not that - 
important . 

Effect of Eesidence Time* The theoretical plot shows 
a gradual decline of M with increasing ~'C , over the range 
of 0.1 hr. to 0.5 hr. as shown in fig. 5. This is to be 
expected since after a certain value of XL is reached, the 
particles become large and due to their increased surface- 
area, absorb more radicals. The increased supply of initiator 
radicals causes small chain to form within the particles and 
thes M n is lowered. Initially of course is '-.zero when C is 
r.ero and must rises to a maximum. 

Then desorption is opera tive, the plot depicts a 
maxima clearly whereafter the curve follows the curve for 
no desorption. Desorption causes size reduction of particles, 
hence the initiator supply to the pai tides is reduced. 
Therefore M increases with C for a longer period than when 
desorption is absent. Once however, the particles have grown 
appreciably, they start absorbing initiator redicals faster 
and the same behaviour of with % is observed as with no 
desorption. So the effect of desorption is to shift the 
maxima towards higher values of £ . 

Fit with the experimental data is poor, but is better 
than that of EeGraff’s theory. The inflexion observed in the 
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Pit of the theoretical curve with experimental data 
is fair, "being about as good as DeGraff's predictions. 

Molecular Weight Distribution : 

The experimental values obtained in the present work 
show a very poor fit with the theoretical prediction of M.Yt.I)., 
due to the faet that the confidence in the experimental results 


is rather low 
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. CHAPTER 7 

f- *. 

CONCLUSIONS 

1 t The present work is an extension of the Thompson and 
Stevens model fis-] for the calculation of number of 
particles, rate of polymerisation and the molecular 

weight distribution, alongwith the absolute particle 
size distribution, , 

2. Results obtained by using this model have been found 
to be similar to those of the RTL model of DeGraff 
Inssome cases the fit of the theoretical prediction with 
experimental values has been found to be better than 
DeGraff's theoretical predictions. 

3. The model could be made better by considering the re- 
entry of live radicals into the polymer particles. 

4. Improvements in the present apparatus would be needed to 
obtain reliable experimental data. In particular proper 
flow rate control is required. 
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Append ix I 


Oheck on the consistency of Boundary Conditions 


We know that: . , 

2/3 


X 


A^y / +■ A^ 


-1/3,. *.*2/3 ) 


, N , % (A 6 y ,/ ' / +A 4 y— 

or f(0,v) = f(l ,v) 1 73 


A A y "+ A c 


Erom (6); at y = 1 , f(l ,v ) = 


f(0,v m ) = f(l ,yJ 


m ' 

a 6 + A 4 
m/ a 4 + a 5 


B i a 11 a m R A 
K 


(la) 


Next, writing a balance on f(o,v) particles lying in a 
volume range of v m to v+ £ v m we have: 


nv 


0 = f (°- 7 m + "2 a > fc v m + f(0 - v m + > &v m a p E ia N A~ f(1 ’ v 


m 




m 


) ^,7 ( Vp + R. a N. \ 
' & m v la p.A) 


Taking limit 4 v 


m 


0 and a. 


a We have:- 

m » 


, s v K o a m , B ia a m N A 

f(o,v m ) = + 


in' v 


m 


1 R- a J* 

— + la m A 


Ar+ A/ , v 

•r f(o,v m ) = f(i,v m > frrl which is same as Ua) 

Hence the consistency' of b.c» on f(l,v m ) has been verified 
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